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Abstract 

The mixed spin- (1/2, 1) Ising chain with axial and rhombic zero- field splitting pa- 
rameters in a presence of the longitudinal magnetic field is exactly solved within 
the framework of decoration-iteration transformation and transfer-matrix method. 
Our particular emphasis is laid on an investigation of the influence of the rhom- 
bic term, which is responsible for an onset of quantum entanglement between two 
magnetic states = ±1 of the spin-1 atoms. It is shown that the rhombic term 
gradually destroys a classical ferrimagnetic order in the ground state and simul- 
taneously causes diversity in magnetization curves including intermediate plateau 
regions, regions with a continuous change in the magnetization as well as several 
unusual field-induced transitions accompanied with magnetization jumps. Another 
interesting findings concern with an appearance of the round minimum in the tem- 
perature dependence of susceptibility times temperature data, the double-peak zero- 
field specific heat curves and the enhanced magnetocaloric effect. The temperature 
dependence of the specific heat with three separate maxima may also be detected 
when driving the system through the axial and rhombic zero-field splitting param- 
eters close enough to a phase boundary between the ferrimagnetic and disordered 
states and applying sufficiently small longitudinal magnetic field. 
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1 Introduction 



Over the last few decades, exactly solved one-dimensional quantum spin mod- 
els [1,2,3] have attracted considerable research interest as they may describe 
subtle quantum phenomena to emerge in real magnetic materials without a 
danger of over-interpretation, which is inherent to any approximative treat- 
ment. The present work is devoted to an exact study of the mixed spin- (1/2, 
1) Ising chain model, which accounts both for the axial zero-field splitting 
(AZFS) as well as the rhombic zero-field splitting (RZFS) parameter in a 
presence of the applied longitudinal magnetic field. It is worthy of notice that 
the special limiting case of this model system in an absence of the external 
magnetic field has been proposed and exactly solved by Wu et al. [4,5,6] using 
the rigorous procedure based on the Jordan- Wigner transformation [7] (see 
Refs. [8,9,10] for related works on this subject). However, it has been recently 
shown by the present authors [11] that the exact results obtained by Wu et 
al. [4,5] can also be recovered by another independent way by making use of 
the transfer-matrix method. The foremost advantage of the formulation based 
on the transfer-matrix method lies in the fact that this rigorous method may 
be even applied in a presence of the non-zero longitudinal magnetic field. The 
main purpose of this work is therefore to investigate the effect of longitudinal 
field on magnetic properties of the mixed spin-(l/2, 1) Ising chain with both 
AZFS and RZFS parameters. 

Before proceeding to an exact calculation for the investigated model system, let 
us briefly comment on an experimental motivation of our study. It is notewor- 
thy that there exist several heterometallic molecular-based compounds with 
a magnetic structure, which can be properly described as one-dimensional 
chain of alternating spin-1/2 and spin-1 metal ions. Among the most common 
examples of the one-dimensional mixed spin-(l/2, 1) chains one could mention 

• CuNi(EDTA).6H 2 0, [12] 

• CuNi(pbaOH)(H 2 0) 3 .nH 2 0, [13] 

• CuNi(pba)(D 2 0) 3 .2D 2 0, [14] 

• PPh 4 [Ni(pn) 2 ][Fe(CN) 6 ]H 2 0, [15] 

• {Pr(bet) 2 (H 2 0) 3 Fe(CN) 6 }, [16] 

• {Ru(acac) 2 (CN) 2 }{Ni(dmphen)(N0 3 )}. [17] 

Even though the vast majority of aforementioned polymeric compounds should 
be preferentially regarded as experimental representatives of the mixed-spin 
quantum Heisenberg chain rather than the mixed-spin Ising chain, it is the 
authors' hope that our exact analytical results for the mixed spin-(l/2, 1) 
Ising chain may provide a deeper insight into several important vestiges of 
real physical behavior and explain some experimental results at least at a 
qualitative level. Besides, one should also expect that the theoretical descrip- 
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tion based on the mixed spin- (1/2, 1) Ising chain may be quite appropriate for 
those heterometallic coordination polymers, where at least one from both con- 
stituent magnetic ions possesses a rather high magnetic anisotropy. It should 
be stressed that the magnetic behavior of this type has been recently found 
in two heterometallic complexes containing highly anisotropic rare-earth ions 
[18,19,20]. 

The outline of the present paper is as follows. In the following Section 2, we 
will shed light on the basic steps of our exact calculation for the investigated 
model system. The Section 3 deals with the discussion of the most interesting 
results obtained for the phase diagrams and basic thermodynamic quantities. 
Finally, some concluding remarks are given in the Section 4. 



2 Exact solution of the mixed-spin Ising chain 



Let us consider the mixed spin-(l/2, 1) Ising chain with AZFS and RZFS 
parameters in a presence of the longitudinal external magnetic field. Suppose 
that the linear chain consist of the alternating spin- 1/2 and spin-1 atoms, 
whereas the former spin-1/2 atoms constitute the sublattice A and the latter 
spin-1 atoms form the sublattice B. The total Hamiltonian of the system can 
be written as a sum of three parts 



Ti — Hex + T~tif s + ^zeei (1) 

which account for the nearest-neighbor Ising interaction, AZFS and RZFS 
terms acting on the spin-1 atoms and the magnetostatic (Zeeman's) energy of 
the spin-1/2 and spin-1 atoms in the applied longitudinal magnetic field 



N 

n ex = -JJ2sm + ^ + i), (2) 

k=l 

N N 



K% = ~D E(^) 2 - E £[(^) 2 - (St)'}, (3) 

k=l k=l 
N N 

H zee = -H A Y J ^~H B Y J S z k . (4) 

k=l k=l 

Above, o z k and S% (a = x, y, z) denote standard spatial components of the 
spin-1/2 and spin-1 operators, respectively, iV denotes a total number of spins 
from each sublattice and the periodic boundary condition <Jn+i = °"i is im- 
posed for simplicity. The parameter J stands for the Ising interaction between 
nearest-neighboring spin-1/2 and spin-1 atoms, whereas the terms D and E 
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label the AZFS and RZFS parameters acting on the spin-1 atoms only. Last, 
two Zeeman's terms Ha and Hb describe the influence of longitudinal mag- 
netic field on the spin- 1/2 and spin-1 atoms from the sublattice A and B, 
respectively. 

Before proceeding further, it is worthwhile to remark that there exist several 
equivalent representations of the zero-field splitting Hamiltonian i-L^fs given by 
Eq. (3). As a matter of fact, one may for instance prove one-to-one correspon- 
dence between %^j s and the effective spin Hamiltonian with three different 
single-ion anisotropy parameters D x , D y and D z 



N N N 

n% = -d* Y,(szf - dv E(^) 2 - d z E(^ 2 ) 2 - (5) 

k=l k=l k=l 

The Hamiltonians Ti^fs an d ^*/s differ one from the other just by unimportant 
constant term, because the mapping relations D = D z — r)X 1 2 DV and E = DX ~ DV 
establish a precise equivalence between the Hamiltonians (3) and (5) (see Ref. 
[11] for more details). It should be also noticed that the particular case of 
the Hamiltonian T-L z f s with D v = has been considered by Wu et al. [4] 
in their recent work. However, it has been shown in our preliminary report 
[11] that the Hamiltonian with one less free parameter is much more 

appropriate for the interpretation of obtained exact results compared with 

~ (2) 

the Hamiltonian 7i z J s and thus, this more convenient definition of the zero- 
field-splitting Hamiltonian will be used throughout the rest of this paper. 

Now, let us turn our attention to the main points of the method, which enables 
an exact treatment of the investigated quantum spin chain. First, the total 
Hamiltonian (1) can be rewritten as a sum of the Zeeman's term for all spin- 
1/2 atoms from the sublattice A and the sum of site Hamiltonians 



N N 

n = J2n k -H A J2° z k , (6) 

k=l k=l 

whereas each site Hamiltonian i-i k involves all the interaction terms including 
the kth. spin-1 atom from the sublattice B 



U k = -E k S z k - D{S z k f - E[{StY - (St) 2 } (7) 

with E k = J(cr k + o'k+i) + B ecause the Hamiltonians (7) at different sites 
commute, i.e. [7-^,7-^] = is valid for each i ^ j, the partition function can 
be partially factorized and consequently rewritten into the form 
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z 



N 

n n ex p 

{<T k }k=l 



'0Hj 



(J, 



Tr 5fc exp(-pH k ) 



N 

En4 

{a fe }fc=l 



(8) 



where /3 = 1/(A; B T), /c B is Boltzmann's constant, T is the absolute tempera- 
ture, Trs k means a trace over degrees of freedom of the fcth spin-1 atom from 
the sublattice B and the symbol J2{v k } denotes a summation over all possible 
configurations of the spin-1/2 atoms from the sublattice A. The crucial step of 
our exact procedure represents calculation of the expression Tr,s fc exp(— fH-i k \ 
For this purpose, it is useful to rewrite the site Hamiltonian (7) into the usual 
matrix representation 



V 



-E k — DO -E 

—E E k -D 



(9) 



using the standard basis of ket vectors \kj) = |±1) , |0) (j = 1 — 3) correspond- 
ing to the three possible spin states S k = ±1,0 of the kth spin-1 atom from 
the sublattice B. The straightforward diagonalization of the site Hamiltonian 
yields the following eigenenergies and eigenvectors 



A W = —D — y[E* + E\ = cos (^) |1> + sin |-1) , 

A fe2 = o, hM = |o), 

\ k3 = -D + y/E 2 k + E*, |^ 3 )=sin(^)|l)+cos(^)|-l), (10) 

with the mixing angle ip k defined through the relation ip k = arctan(-J^). It 
is worth mentioning that the eigenenergies listed in the set of Eqs. (10) can 
readily be used for calculating the expression Tr 5fe exp(— (31-ik) and moreover, 
the analytical form of the site partition function Z k then immediately implies 
a possibility of performing the generalized decoration-iteration transformation 
[21,22,23] 



Z k = exp 
= A exp 



+ <7fc+i) 



1 + 2 exp(/3D) cosh py/E% + E 2 



(ii) 



The unknown mapping parameters A, R and H entering in Eq. (11) can be 
obtained from the self-consistency condition of the decoration-iteration trans- 
formation. Following the standard procedure [21,22,23], one directly obtains 
explicit expressions for all three transformation parameters 
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A=(V 1 V 2 V 3 2 )*, 



f3H = In 




(12) 



which are defined in terms of newly defined functions Vi, V2 and V3 




expl^j 1 + 2 exp(/3D) cosh (J + H B ) 2 + E 2 



1 + 2 exp(/3D) cosh (f3^J(J - H B ) 2 + E 2 



V 3 



1 + 2 exp(/3D) cosh U^JhI + E 2 



(13) 



The functions Vi, V2 and V3 physically correspond to three independent ex- 
pressions for the site partition function Z k , which can be obtained from the 
transformation formula (11) by considering all four spin configurations of two 
Ising spins a k and cr k+1 involved therein (recall that they enter into the param- 
eter E k as well). The substitution of the decoration-iteration transformation 
(11) into the formula (8) leads in turn to a rigorous mapping relationship be- 
tween the partition function Z of the mixed-spin Ising chain and the partition 
function Z ic of the simple spin-1/2 Ising chain with the effective temperature- 
dependent nearest-neighbor interaction R and magnetic field H 



The relationship (14) between both partition functions proves that there exist 
a rigorous mapping correspondence between the mixed-spin Ising chain and 
the simple spin-1/2 Ising linear chain, since both partition functions differ one 
from another just by the multiplicative factor A given by Eqs. (12)-(13). In 
this regard, the equality (14) may be considered as a mathematical proof of the 
aforementioned mapping equivalence and hence, the partition function of the 
mixed-spin Ising chain can exactly be calculated from the well-known formula 
for the partition function of the simple spin-1/2 Ising linear chain [1] with the 
effective nearest-neighbor interaction R and magnetic field H unambiguously 
determined by Eqs. (12)-(13). Besides, the mapping relationship (14) is also 
very suitable for deriving exact results for other quantities, such as free energy, 
magnetization, entropy, specific heat, susceptibility, etc. For instance, the free 
energy of the mixed-spin Ising chain per unit cell reads 



At this stage, both sublattice magnetizations can easily be calculated as a 
partial derivative of the the free energy with respect to the relevant magnetic 



Z(f3, J, D, E, H A , H B ) = A N Z lc (f3, R, H ). 



(14) 




(15) 




After straightforward, 
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but a little bit cumbersome arrangement, both sublattice magnetizations can 
be written in this compact form 



1 



m A = (ZD = 77 



v 1 -v 2 



2 ^J(V 1 -V 2 ) 2 + W 2 



(16) 



m B = (S z k 



(vi - v 2 ){w 1 - w 2 ) + 4p 3 + (w 1 + w 2 y04 - v 2 y + wi 



(v, - v 2 y + w 2 + (v, + v 2 ) Jiy 1 - v 2 



,(17) 



where (• • •) means the standard canonical ensemble averaging and the func- 
tions Wi, W 2 and W3 are defined as follows 



Wi = 2 
W 2 = -2 



( J + ^ :expU + ^smh 



J + H B ) 2 + E 2 



^{J + H B ) 2 + E 2 



{J-H_ 



B 



) 2 + E 2 



exp ( (3D — — — —J sinh 



exp(/3D) sinh pJH B + E 2 



(3^{J-H B ) 2 + E 2 
(18) 



The total magnetization normalized with respect to its saturation value is 
then given by m/m s = 2(mi + mB)/3. Other basic thermodynamic quantities 
such as entropy and specific heat can also be readily calculated from the free 
energy (15) with the help of basic thermodynamic relations 



S = 



dT 



C = -T 



11 



'd 2 F s 

dr 2 



(19) 



H 



Notice that the final expressions for the entropy and specific heat obtained by 
the use of Eq. (19) are too cumbersome to write them down here explicitly. 
Furthermore, one may also obtain the relevant exact result for the initial 
longitudinal susceptibility reduced per elementary unit from the previously 
derived exact results (16)— (17) for both sublattice magnetizations 



Xz 



dm/ 



+ 



'dm 



B 



T,H A ^0 



0H f 



= 



Ui , ui + u 2 (u 4 + u 5 + u 6 ) 



+ 



c/ 2 (c/i + c/ 2 ) 



(20) 



For simplicity, the functions [/& (k — 1 — 6) to emerge in Eq. (20) for the initial 
longitudinal susceptibility are explicitly given in the Appendix. 
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On the other hand, the exact results derived for the free energy (15) or both 
sublattice magnetizations (16)— (17) cannot be utilized for calculating the ini- 
tial transverse susceptibility. For this purpose, one first needs to calculate the 
relevant spatial component of the quadrupolar moment in the limit of zero 
external magnetic field (Ha = H B = 0.0) and to relate it subsequently to the 
initial transverse susceptibility using the fluctuation-dissipation theorem [24]. 
It is noteworthy that all three spatial components of the quadrupolar moment 
can easily be obtained with the help of the generalized Callen-Suzuki identity 
[25,26]. According to this, the statistical mean values determining the spatial 
components of the quadrupolar moment of the spin-1 atoms may be calculated 
from the exact spin identity 



a\2 



TT Sk [(St?eM-m k )] 
Tr 5fe exp(-/374) 



[a 



x,y,z) 



(21) 



After straightforward but rather lengthly calculation based on the exact spin 
identity (21) and the relevant eigenvalues and eigenvectors (10) of the site 
Hamiltonian (7), all three spatial components of the quadrupolar moment can 
be expressed through the following formulas 



<f B 



ut + u} 

U, + U 2 ' 



Qb 



Ui + u 2 
C/i + u 2 ' 



z U x + U 2 - 2 
Qb = 



U ± + U 2 



(22) 



where the functions Ui and U\ are defined for i — 1, 2 and j — ± in Appendix. 
Now, both spatial components of the initial transverse susceptibility reduced 
per elementary unit can be easily calculated by making use of the fluctuation- 
dissipation theorem 




■ N 

k=i 



i 2, 



■ N 

E(3? + *Z) 
.k=i 



in = x,y) (23) 



In the zero-field limit, the exact expression for the initial transverse suscep- 
tibility resulting from the fluctuation-dissipation theorem largely simplifies, 
because all nearest-neighbor as well as further-neighbor pair correlation func- 
tions vanish in an absence of pair interactions between transverse (x or y) 
spatial components of the spin operators. Both spatial components of the 
initial transverse susceptibility can be thus expressed solely in terms of the 
respective spatial components of the quadrupolar moment 



X 7 = P 



1 



+ ((^ 7 ) 2 







(24) 
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which are exactly known from Eq. (22) derived previously using the generalized 
Callen-Suzuki identity [25,26]. 



3 Results and discussion 



In this section, we will focus our attention to the most interesting results to be 
obtained for phase diagrams and basic thermodynamic quantities. It is worth- 
while to remark that the magnetic behavior of the investigated model system 
in an absence of the external magnetic field has been thoroughly examined 
in Refs. [4,5,11] and thus, the respective behavior in a non-zero longitudinal 
magnetic field will be the main focus of the present paper. It should be also 
noticed that all results obtained in the preceding section are quite general as 
they hold both for the special case with the ferromagnetic (J > 0) as well as 
antiferromagnetic (J < 0) interaction. The ferromagnetic interaction between 
the nearest-neighbor spins generally leads to the ferromagnetic spin alignment, 
while the antiferromagnetic interaction causes the ferrimagnetic spin arrange- 
ment. It could be expected that the magnetic behavior of the ferrimagnetic 
model in the applied longitudinal magnetic field should be much more diverse 
compared with its ferromagnetic counterpart and therefore, we will henceforth 
assume J < for the sake of simplicity. In order to reduce the total number 
of free Hamiltonian parameters, our subsequent analysis will be restricted to 
the special case with the uniform magnetic field Ha = Hb = H acting on the 
spin-1/2 and spin-1 atoms, which physically corresponds to the situation with 
the equal g-factors for both kinds of the magnetic atoms. 



3.1 Ground- state phase diagram 

Let us construct first the ground-state phase diagram, which will clarify the 
magnetic behavior of the ferrimagnetic mixed spin- (1/2,1) Ising chain with the 
AZFS and RZFS parameters in the longitudinal magnetic field. The ground- 
state phase diagram in the D / \ J\— Hj \ J\ plane is illustrated in Fig. 1(a) for 
two different values of the RZFS parameter Ej \ J\ = 0.0 and 0.5. As one can 
see, the investigated mixed-spin chain keeps under the influence of longitudinal 
magnetic field one from three available ground states, which are separated one 
from another by discontinuous (first-order) phase transition lines. Three spin 
arrangements emerging within the available ground states can be described 
through the following eigenvectors. 
Quantum ferrimagnetic phase: 
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(a) D/\J\ (b) D/IJI 

Fig. 1. The ground-state phase diagram in the D / \ J\ —H/ \J\ plane of the ferrimag- 
netic (J < 0) mixed spin-(l/2,l) Ising chain for: (a) two different values of RZFS 
parameter Ej \J\ = 0.0 (solid lines) and Ej \J\ = 0.5 (solid lines with symbols), (b) 
the detail of the particular case with Ej \J\ = 0.5 in the vicinity of the triple point. 



N 



\fp) = n 

k=l 



er t = — 



ip± = arctan 



1 

2 

E 



cos 



fx 
2 



\S* = 1) + sin 



2 



I Q z 



■i) 



H 



B 



J 



(25) 



Disordered phase: 



\dp) = n 

k=l 



o2 = 8gn(f0-)|52 = 0>, 



(26) 



Quantum paramagnetic phase: 



N 



\pp) = n 

k=l 



■ if) 1^ = 1)+ sin (f )\S' k = -l) 



(^2 = arctan 



E 



H 



B 



\J\ 



(27) 



The spin alignment becomes quite obvious when assuming the particular case 
without the RZFS anisotropy, i.e. Ej \ J\ = 0.0. The eigenvector \FP) then 
refers to the classical ferrimagnetic phase with a perfect antiparallel orientation 
of the nearest-neighboring spin-1/2 and spin-1 atoms each of them occupy- 
ing the spin state |cr| = —1/2) and |S| = 1), respectively. This phase appears 
in the ground state just if the AZFS parameter is greater than the bound- 
ary value D b / \ J\ = —1.0, otherwise the disordered phase \DP) becomes the 
ground state. In an absence of the external magnetic field, each spin-1/2 atom 
randomly occupies in \DP) any of its two available spin states \o z k = ±1/2) on 
behalf of the non-magnetic character |S| = 0) of all the spin-1 atoms, which 
is enforced by the AZFS parameter. In any non-zero field, all the spins- 1/2 
atoms tend to align into the external-field direction as they effectively behave 
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as free paramagnetic spins within \DP). Finally, the classical paramagnetic 
phase \PP) becomes the ground state whenever the external field exceeds the 
saturation value HJ \J\ = 2.0 for D/ \J\ > -1.0 or HJ \ J\ = 1 - D/\J\ for 
D/\J\ < —1.0. The classical paramagnetic phase can be characterized by a 
perfect alignment of all the spin-1/2 as well as spin-1 atoms into the external- 
field direction. All three possible ground states coexist together at a triple 
point given by the coordinates \D t j \J\ ,H t / |J|] = [—1.0,2.0]. 

Contrary to this, the situation becomes much more involved when assuming 
the non-zero RZFS parameter and its influence on the ground-state spin ar- 
rangement. The typical ground-state phase diagram is shown in Fig. 1(a) for 
one specific value of the RZFS parameter E/ | J| = 0.5. The non-zero RZFS 
term introduces into the otherwise classical mixed-spin Ising chain x- and y- 
components of the spin operators and hence, this parameter is responsible for 
the onset of local quantum fluctuations that consequently lead to a quantum 
entanglement between two magnetic spin states \S%. = ±1) of the spin-1 atoms. 
In this respect, the RZFS parameter gradually destroys the perfect ferrimag- 
netic spin arrangement, which exists in \FP) in an absence of the RZFS term. 
Indeed, the quantum superposition of the spin states \S% = 1} and \S% = — 1) 
governs the overall magnetic behavior of the spin-1 atoms in \FP), whereas 
the relevant spin state |<r| = —1/2) of the spin-1/2 atoms is not affected by 
the RZFS term at all. It is also quite evident from the mixing angle tpi given 
by Eq. (25) that the RZFS term generally enhances the probability amplitude 
of the less probable spin state |S| = —1) on account of the more probable spin 
state |<Sfc = 1), while the reverse is the case when considering the effect of the 
longitudinal magnetic field. If the external field exceeds some critical value, the 
investigated spin system passes from the quantum ferrimagnetic phase \FP) 
to the quantum paramagnetic phase \PP). All the spin-1/2 atoms tend to 
align towards the external- field direction in \PP), while the quantum entan- 
glement of the spin states \S% — 1) and \S% = —1) persists in this phase even 
if the rising magnetic field gradually suppresses the quantum superposition 
defined through the another mixing angle ip2 given by Eq. (27). 

As far as the disordered phase \DP) is concerned, there is no change in the 
relevant spin arrangement of this phase in comparison with the previously 
discussed specific case without the RZFS term. Namely, the RZFS parameter 
does not affect the overall spin arrangement inherent to \DP), since this term 
is responsible only for the quantum entanglement between two magnetic spin 
states \Sl = ±1), whereas all the spin-1 atoms reside the non-magnetic spin 
state \S% — 0) in \DP). However, there is an interesting shift of the ground- 
state phase boundary between \FP) and \DP) to a more negative values 
of the AZFS parameter, which emerges as a result of the influence of the 
RZFS parameter. In fact, this phase boundary becomes a striking non-linear 
curve instead of a simple vertical line for any non-zero RZFS parameter [see 
the detail depicted in Fig. 1(b) as well]. In a relatively narrow range of the 
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Fig. 2. Both sublattice magnetizations as well as the total magnetization re- 
duced with respect to its saturation value versus the external magnetic field 
at the sufficiently low temperature ksT/\J\ = 0.01, the fixed value of the 
RZFS parameter E/\J\ = 0.5 and four different values of the AZFS parameter 
D/\J\ = -0.6,-1.4,-1.07,-1.1. 

AZFS parameters, e.g. D / \ J\ G (—1.115, —1.040) for the particular case with 
E 1 1 J\ = 0.5, the investigated mixed-spin chain thus displays a peculiar se- 
quence of two field-induced phase transitions. At lower critical field, the system 
first undergoes the field- induced transition from \FP) to \DP), while another 
field- induced transition from \DP) to \PP) takes place at upper critical field. 
The triple point, at which all three phases coexist together in the ground state, 
is given by the coordinates [D t / \ J\ ,H t / |J|] = [—1.040,1.899] for the specific 
case with E / \ J\ = 0.5. 



3.2 Magnetization process 



Now, let us explore in detail the effect of AZFS and RZFS parameters on the 
magnetization process. Both sublattice magnetizations (m^, mg) as well as 
the total magnetization reduced with respect to its saturation value {m/m s ) 
are plotted in Fig. 2(a)-(d) against the external magnetic field for one se- 
lected value of the RZFS parameter Ej \ J\ = 0.5 and four different values of 
the AZFS parameter. The value D / \ J\ = —0.6 is chosen so as to achieve the 
quantum ferrimagnetic (paramagnetic) phase \FP) (\PP)) in the ground state 
at low (high) enough magnetic fields. In this particular case, the total magneti- 
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zation at first rises steadily as the longitudinal magnetic field strengthens until 
it reaches a certain critical field at which one observes the magnetization jump 
due to the abrupt spin reversal of all the spin- 1/2 atoms into the external- 
field direction. The initial continuous increase of the total magnetization as 
well as the one observed after the critical field can evidently be attributed 
to a gradual increase of the occurrence probability of the majority spin state 
\Sj: = 1), which appears in \FP) or \PP) at the expense of the occurrence 
probability of the minority spin stcttc \Sfc — — 1} owing to the magnetic field 
strengthening. It is noteworthy that the strict magnetization jump serving in 
evidence of the field- induced transition from \FP) to \PP) takes place just at 
zero temperature, however, the rather steep (but still continuous) increase of 
the magnetization in a close neighborhood of the critical field reflects this zero- 
temperature transition at sufficiently low (but non-zero) temperatures as well. 
Furthermore, the second selected value of the AZFS parameter D / \ J\ = —1.4 
is chosen strong enough in order to cause an appearance of the disordered 
(paramagnetic) phase \DP) (\PP)) in the ground state at sufficiently low 
(high) magnetic fields. In this special case, the total magnetization exhibits 
at low magnetic fields the intermediate magnetization plateau at one third of 
the saturation magnetization, which only comes from a perfect alignment of 
the spin-1/2 atoms into the external-field direction (the spin-1 atoms do not 
contribute to the total magnetization at all as they reside the non-magnetic 
state \S%. = 0) within \DP)). At a certain critical field, the quantum spin chain 
undergoes the field-induced transition from \DP) to \PP), which is accompa- 
nied with the abrupt change in the spin state of the spin-1 atoms from the 
non-magnetic state \Sj: = 0) to the entangled pair of states \Sj: = ±1). In the 
consequence of that, the total magnetization rises steadily with the magnetic 
field after the relevant critical field owing to a gradual change in the occurrence 
probabilities of the majority and minority spin states \S%. = 1) and \S% = — 1), 
which appears in \PP) under the influence of the longitudinal magnetic field. 

It should be nevertheless mentioned that the most interesting magnetization 
curves can be detected if a relative strength of the AZFS parameter is se- 
lected from the relatively narrow interval D/\J\ G (—1.115,-1.040) when 
assuming Ej \ J\ = 0.5. Figs. 2(c), (d) provide two illustrative examples of the 
low-temperature magnetization curves with the AZFS parameter taken from 
this range, which ensures occurrence of two sequential field-induced transitions 
from \FP) to \DP) and subsequently from \DP) to \PP) upon the magnetic 
field strengthening. Actually, it can be clearly seen from Figs. 2(c),(d) that the 
total magnetization at first continuously increases within \FP) with increasing 
the magnetic field until it reaches the lower critical field, where the magneti- 
zation jumps towards the intermediate plateau with the constant value of the 
total magnetization that is stable for mediate values of the external magnetic 
field. The plateau region, which corresponds to the appearance of \DP), then 
ends up at the upper critical field reflecting another field-induced transition 
from \DP) to \PP). Accordingly, the total magnetization jumps from one third 
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Fig. 3. The total magnetization reduced with respect to its saturation value as a 
function of the magnetic field at the sufficiently low temperature ksT/ \J\ = 0.01, 
one specific value of the RZFS parameter E/\J\ = 0.5 and several values of the 
AZFS parameter: (a) D/\J\ = -0.6, -1.2, -1.4; (b) D/ \J\ = -1.05, -1.07, -1.1. 

of the saturation value to a rather high value relatively close to the satura- 
tion magnetization and then, it rises steadily with increasing the longitudinal 
magnetic field within \PP) until it reaches the saturation magnetization in 
the limit of infinitely strong magnetic field H/ \J\ — > oo. 



For the sake of comparison, the total magnetization reduced with respect to its 
saturation value is plotted against the external magnetic field in Figs. 3(a),(b) 
for one specific value of the RZFS parameter Ej \ J\ = 0.5 and several values 
of the AZFS parameter at low enough temperature k^T '/ ' \J\ = 0.01. These 
figures clearly demonstrate a possible diversity in the magnetization process, 
which might include the regions with the continuous change in the total mag- 
netization, the intermediate magnetization plateau, as well as, one or two 
magnetization jumps accompanying the field-induced transitions between dif- 
ferent phases. Finally, it should be also noted here that the temperature rise 
gradually smoothens the magnetization behavior close to and at the intermedi- 
ate magnetization plateau, which cannot be consequently seen in the relevant 
magnetization curves at sufficiently high temperatures. 



3.3 Specific heat 



Another quantity, which is important for overall understanding of thermody- 
namics, is the magnetic contribution to the specific heat. Temperature varia- 
tions of the zero- field specific heat are shown in Figs. 4(a), (b) for the RZFS 
parameter Ej \J\ = 0.5 and several values of the AZFS parameter. As one 
can see, the temperature dependence of the specific heat always exhibits at 
least one round Schottky-type maximum regardless of whether \FP) or \DP) 
constitutes the ground state. The specific heat has a single rounded maxi- 
mum when considering either positive or small negative values of the AZFS 
parameter warranting the appearance of \FP) in the ground state such as 
D/ | J | = 0.0 and —0.6. If the AZFS parameter is selected close enough to the 
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Fig. 4. Temperature changes of the specific heat for one selected value of the RZFS 
parameter E/\J\ = 0.5 and several values of the AZFS parameter. Figs. 4(a)-(b) 
show thermal variations of the specific heat in the zero-field case H/\J\ = 0.0, 
whereas the special case D / \ J\ = — \/l.25 corresponds to the phase boundary be- 
tween \FP) and \DP). Figs. 4(c)- (d) illustrate thermal variations of the specific 
heat in the non-zero magnetic field when assuming D/ \ J\ = —1.3. 

phase boundary between \FP) and \DP), then, the additional second peak 
can be observed at low temperatures no matter whether \FP) or \DP) con- 
stitutes the ground state (see the specific heat curves for D / \ J\ = —1.05 and 
—1.2). However, the pronounced double-peak dependence of the specific heat 
gradually changes into the dependence with a single more or less symmetric 
maximum when considering more negative values of the AZFS parameter far 
from the phase boundary between \FP) and \DP) (see the thermal depen- 
dences for D 1 1 J\ = —1.4 and —1.8). These observations would suggest that 
the low-temperature peak in the double-peak dependence of the zero-field 
specific heat arises from thermal excitations between the ground-state spin 
configuration and the relevant excited state with only slightly higher energy. 

The situation becomes even more intriguing when assuming the AZFS param- 
eter close enough to the phase boundary between \FP) and \DP) and applying 
the non-zero magnetic field. It is quite apparent from Fig. 4(c) that the ther- 
mal dependence of specific heat with a single more or less symmetric maximum 
can be recovered upon strengthening the external magnetic field (see the curve 
for H/ \ J\ = 0.5). However, one may also detect the very unusual triple-peak 
dependence of the heat capacity when applying the sufficiently small but non- 
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(a) k B T/\J\ (b) k B T/\J\ 



Fig. 5. The temperature dependence of the initial longitudinal susceptibility 
times temperature (XzT) product for one selected value of the RZFS parame- 
ter E/\J\ = 0.5 and several values of the AZFS parameter. The particular case 
D/\J\ = — \/l-25 corresponds to the phase boundary between \FP) and \DP). 

zero longitudinal magnetic field as depicted in Fig. 4(d). It is noteworthy that 
the third peak is located in the low-temperature area and the emergence of 
this additional third peak can be explained as the Zeeman's splitting of energy 
levels. In agreement with this argument, the third peak shifts linearly towards 
higher temperatures with increasing the strength of the magnetic field until it 
completely merges with the second low-temperature maximum. Therefore, the 
triple-peak dependence of the heat capacity may be observed in a relatively 
weak magnetic fields only. 

3.4 Longitudinal susceptibility 

Next, let us examine the thermal dependence of the initial longitudinal sus- 
ceptibility times temperature product XzT as shown in Fig. 5. If the AZFS 
parameter is selected so as to achieve \FP) in the ground state, then, \ Z T 
product generally exhibits an exponential divergence when approaching zero 
temperature as depicted in Fig. 5(a). It can be seen from this figure that \ Z T 
monotonically decreases with increasing temperature for positive or small neg- 
ative values of the AZFS parameter. In addition, one may also detect a more 
striking non-monotonous temperature dependence of XzT with a rather flat 
minimum when the AZFS parameter is chosen sufficiently close but slightly 
above the ground-state boundary between \FP) and \DP) [see for instance 
the curve D/\J\ = —1.1 in Fig. 5(a)]. The emergence of the round mini- 
mum in the temperature dependence of x z T product is a typical feature of 
the quantum ferrimagnets [27], because the monotonous decrease (increase) of 
X Z T product with increasing temperature indicates predominantly ferromag- 
netic (antiferromagnetic) character of excitations. The position of the round 
minimum can be thus regarded as a crossover point from the region with pre- 
dominant ferromagnetic towards the region with predominant antiferromag- 
netic excitations. If the AZFS parameter is strong enough to stabilize \DP) 
in the ground state [see Fig. 5(b)], then, \ Z T product increases monotonically 
with increasing temperature from its minimum initial value 1/4. This constant 
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value evidently comes from the response of the spin- 1/2 atoms with respect to 
the infinitesimal change of the magnetic field, because the spin-1 atoms occupy 
within the \DP) the non-magnetic spin state \S% — 0). It is also noteworthy 
that XzT product converges to the constant value of 11/12 in the limit of high 
temperatures T — > oo regardless of whether \FP) or \DP) is being the ground 
state. This value agrees with the Curie law for the mixed spin-(l/2,l) system, 
which consist of an equal number of independent spins-1/2 and spin-1 atoms. 

3.5 Transverse susceptibility 

Further, let us proceed to a discussion of temperature dependences of the 
initial transverse susceptibility times temperature product. It is worthwhile to 
recall that the RZFS parameter effectively introduces a magnetic anisotropy 
into both transverse (x and y) spatial components of the mixed-spin Ising 
chain refined by the AZFS and RZFS parameters and consequently, thermal 
variations of both spatial components of the initial transverse susceptibility 
(Xx an d Xy) might significantly differ one from another. Note furthermore that 
the positive value of the RZFS parameter corresponds to the situation when 
the x-axis is easier magnetization axis than the y-axis within the xy-plane. 

First, let us compare thermal dependences of both initial transverse suscep- 
tibilities by considering the ground-state spin arrangement inherent to \FP) 
(the left panel in Fig. 6). It is noteworthy that both spatial components of 
the initial transverse susceptibility exhibit very different temperature varia- 
tions of the susceptibility times temperature product in this particular case. 
X y T product initially decreases to its local minimum and then gradually in- 
creases with increasing temperature for positive or small negative values of the 
AZFS parameter (e.g. D/\J\ = 1.0 and 0.0), while it monotonously increases 
upon temperature increase for more negative values of the AZFS parame- 
ter (e.g. D/\J\ = —0.6 and —0.9). On the other hand, XxT product always 
exhibits a non-monotonous temperature dependence with a single rounded 
maximum. Another interesting observation is that XxT is in general greater 
than the corresponding value of XyF irrespective of temperature. According 
to Eq. (24), the greater (smaller) value of XxT {XyT) relates to the greater 
(smaller) value of the x(|/)-component of the quadrupolar moment q% (q y B ). 
The zero-temperature limits of the transverse susceptibility times temperature 
product are actually consistent with the asymptotic values of the quadrupolar 
moment obtained for \FP) from Eq. (22) in the zero-temperature limit T — > 



This result is taken to mean that the x-axis is indeed easier magnetization 




(28) 
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Fig. 6. Thermal variations of the initial transverse susceptibility times temperature 
product for the zero longitudinal magnetic field H/\J\ = 0.0, one particular value 
of the RZFS parameter E/\ J\ = 0.5 and several values of the AZFS parameter. The 
upper (lower) panel shows the transverse susceptibility for the x-axis (y-axis), the 
left (right) panel illustrates the dependences for the ground state \FP) (\DP)). 

axis than the y-axis in \FP). It is worthy to notice, moreover, that the zero- 
temperature value of transverse spatial component of the quadrupolar moment 
(28) has been confused with the zero-temperature value of the transverse sus- 
ceptibility in the recent work of Wu et al. [4]. 

It should be also mentioned that the difference in thermal dependences of 
both initial transverse susceptibilities becomes much less pronounced when 
considering the mixed spin- (1/2,1) Ising chain in the another possible ground 
state \DP) (the right panel in Fig. 6). The spin-1 atoms occupy within \DP) 
the non-magnetic spin state \S% — 0), which leads in turn to the equality of 
both transverse spatial components of the quadrupolar moment simultane- 
ously achieving their maximum value q% = q v B = 1.0 in the ground state. 
Owing to this fact, both spatial components of the initial transverse suscepti- 
bility times temperature product start from its maximum value 5/4 and they 
must generally decrease with increasing temperature due to the temperature- 
induced lowering of the quadrupolar moment. The only non- monotonous tem- 
perature dependence of the initial transverse susceptibility times temperature 
product can be found in \ X T vs. T dependence if one considers the AZFS 
parameter close enough to the phase boundary between \DP) and \FP). In 
this special case, \ X T initially decreases rather rapidly to its local minimum, 
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(c) H/\J\ (d) 

Fig. 7. Adiabatic changes of the temperature with the magnetic field for one selected 
value of the RZFS parameter E/\J\ = 0.5 and four different values of the AZFS 
parameter, which correspond to the magnetization curves depicted in Fig. 2. 

then it rises steadily to a shallow local maximum and finally, it repeatedly 
exhibits a gradual decline with increasing temperature. 




3.6 Adiabatic demagnetization 



Last but not least, let us examine the adiabatic demagnetization in connec- 
tion with a possibility of observing the enhanced magnetocaloric effect that 
might be of technological relevance for magnetic refrigeration [28,29,30,31,32]. 
To gain an insight into the adiabatic demagnetization, Fig. 7 displays the 
temperature as a function of the external magnetic field under the adiabatic 
condition (i.e. the constant entropy) for one specific value of the RZFS param- 
eter E/ \ J\ = 0.5 and four different values of the AZFS anisotropy. Note fur- 
thermore that the displayed dependences correspond to the low-temperature 
magnetization curves depicted in Fig. 2. The observed temperature changes 
can easily be understood from a comparison of Figs. 1 and 7. If the AZFS pa- 
rameter is so selected that \FP) constitutes the zero- field ground state [Fig. 
7(a)], then, the spin system finally approaches a relatively low but non-zero 
temperature in the limit of vanishing external magnetic field. Contrary to this, 
the temperature tends towards the absolute zero when a sufficiently strong 
(negative) value of the AZFS parameter energetically stabilizes the disordered 
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(a) H/\J\ (b) H/\J\ 



Fig. 8. Adiabatic changes of the temperature with the external magnetic field for 
one selected value of the RZFS parameter Ej \J\ = 0.5 at relatively low values of 
the entropy. The small dip in the temperature versus the magnetic field dependence 
corresponds to the field-induced transition between \FP) and \DP). 

ground state \DP) in the zero-field limit H/\J\ — > 0.0 [Figs. 7(b)-(d)]. From 
this point of view, the mixed spin-(l/2,l) Ising chain becomes a more effi- 
cient refrigerant if the AZFS parameter drives the system into the disordered 
ground state \DP). It is also worthy to notice that the character of mag- 
netocaloric effect basically changes at transition fields, where one generally 
observes a rather fast cooling (heating) when approaching the relevant tran- 
sition field from above (below). Besides, it is also quite apparent from Fig. 7 
that the relatively fast and efficient cooling (heating) can be achieved only for 
the adiabatic processes with the constant entropy close enough to the value 
S = Nk-Qln2. Under this circumstance, the temperature falls rather quickly 
to zero if the external magnetic field is approaching either the transition field 
between \FP) and \PP), \DP) and \PP), or it tends towards zero. On the 
other hand, the temperature reaches just some non-zero (even if relatively 
small) value if the external magnetic field approaches the transition field be- 
tween \FP) and \DP). It is noteworthy that the latter behavior cannot be 
clearly seen within the displayed scale from Fig. 7 and thus, it is separately 
demonstrated in Fig. 8. It is quite evident from this figure that a relatively 
small dip (minimum) can be observed at the field-induced transition between 
\FP) and \DP), whereas this local minimum in the magnetic field dependence 
of the temperature already disappears upon a small increase of the entropy of 
the adiabatic process. 



4 Concluding remarks 

In this article, the mixed spin-(l/2,l) Ising chain with the AZFS and RZFS 
parameters in the applied longitudinal magnetic field has been exactly solved 
by combining the decoration-iteration transformation with the transfer-matrix 
method. Exact results for phase diagrams, free energy, magnetization, quadrupo- 
lar moment, entropy, specific heat, longitudinal and transverse initial suscep- 
tibilities have been derived and discussed in detail. It has been shown that 
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the RZFS parameter basically modifies the magnetic behavior of the investi- 
gated mixed-spin chain, since it is responsible for the quantum entanglement 
between two magnetic states |<S| = ±1) of the spin-1 atoms. 

Among the most interesting findings to emerge from the present study one 
could mention an extraordinary diverse magnetization process, which may 
include the regions with the continuous change in the total magnetization, 
the intermediate magnetization plateau at one third of the saturation mag- 
netization, as well as, one or two magnetization jumps accompanying the 
field-induced transitions between different phases. Another interesting find- 
ings concern with an appearance of the round minimum in the temperature 
dependence of susceptibility times temperature data, the double-peak zero- 
field specific heat curves and the enhanced magnetocaloric effect in a vicinity 
of the field-induced phase transitions. It has been evidenced that the triple- 
peak temperature dependence of the specific heat may also be found when 
driving the system through the AZFS and RZFS parameters close enough to 
a phase boundary between \FP) and \DP) and applying sufficiently small but 
non-zero longitudinal magnetic field. 



5 Appendix 



The list of functions, which enter into the exact analytical formulas (20) and 
(22) for the initial longitudinal susceptibility and the quadrupolar moment in 
the limit of the vanishing longitudinal magnetic field Ha = Hb — > 0: 



E/i = 1 + 2 exp(l3D) cosh J 2 + E 2 ) , 
U 2 = 1 + 2 exp(/3£>) cosh (pE) , 

U 3 = exp(PD) sinh (pv^PTW) , 

U, = p {j2 2 f E2)3/2 exp(/3D) sinh (fiVJ^+E^) , 
U 5 = J2 + E2 exp(/3D) cosh (j3y/ J 2 + E 2 ) , 
t/ 6 = ^exp(/3D)sinh(/3£), 

Uf = l+ exp{f3D) cosh (fiVJ 2 + E 2 ) ± ^7== s inh (N J 2 + E 2 ) , 
U£ = l + exp(/3D) cosh (f3E) ± exp(/3D) sinh (f3E) . 
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